Local defect in a magnet with long-range interactions 
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We investigate a single defect coupling to the square of the order parameter in a nearly critical 
magnet with long-range spatial interactions of the form r ~^ d+cr \ focusing on magnetic droplets 
nucleated at the defect while the bulk system is in the paramagnetic phase. To determine the static 
droplet profile, we solve a Landau-Ginzburg- Wilson action in saddle point approximation. Because 
of the long-range interaction, the droplet develops a power-law tail which is energetically unfavorable. 
However, as long as a > 0, the tail contribution to the droplet free energy is subleading in the limit 
of large droplets; and the free energy becomes identical to the case of short-range interactions. We 
also consider the effects of fluctuations and find that they do not change the functional form of the 
droplet as long as the bulk system is noncritical. Finally, we study the droplet quantum dynamics 
with and without dissipation; and we discuss the consequences of our results for defects in itinerant 
quantum ferromagnets. 

PACS numbers: 75.10.Jm, 75.10.Nr, 75.40.-s 



I. INTRODUCTION 

In a nearly critical system, a local defect that prefers 
the ordered phase can induce the nucleation of a droplet 
of local order in the nonordered background. Such 
droplets arise, e.g., in disordered systems due to the pres- 
ence of rare strongly coupled spatial regions. They can 
have surprisingly strong consequences for the properties 
of the phase transition. In a classical magnet at nonzero 
temperature, a large finite-size droplet does not have a 
static magnetization; instead it fluctuates very slowly be- 
cause flipping the droplet requires coherently changing 
the order parameter in a large volume. More than 30 
years ago, Griffiths 1 - showed that rare regions and the 
magnetic droplets formed on them, lead to a singularity 
in the free energy in a whole temperature region above 
the critical point, which is now known as the Griffiths re- 
gion or the Griffiths phased Later, it was shown that this 
singularity is only an essential one 3,4,5 and thus probably 
unobservable in experiment (see also Ref. @). 

The effects of magnetic droplets are greatly enhanced if 
the underlying defects are extended macroscopic objects 
(linear or planar defects). In these cases, the droplet dy- 
namics is even slower and so increases their effects. This 
was first found in the McCoy- Wu model^ a 2D Ising 
model with linear defects. Later it was studied in great 
detail in the context of the quantum phase transition of 
the random transverse-field Ising model where the de- 
fects are extended in the imaginary time dimension 
In these systems, the Griffiths singularity in the free en- 
ergy actually takes a power-law form, and the suscepti- 
bility diverges inside the Griffiths region. 

Recently, it has been shown that even stronger effects 
than these power-law Griffiths singularities can occur in 
Ising magnets with planar defects i 11 ! 12 Droplets that are 
extended in two dimensions can undergo the magnetic 
phase transition (and develop a static order parameter) 
independently from the bulk system. This leads to a de- 
struction of the global sharp phase transition by smear- 
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ing. (Note that an unusual magnetization-temperature 
relation was already found in the numerical mean-field 
analysis, Ref. [ID, but it was interpreted as power-law 
critical behavior with a very large exponent /?.) Similar 
smeared phase transitions have also been found in a non- 
equilibrium system in the presence of linear defects 1 ^. A 
recent review of these and other rare region effects can 
be found in Ref. [TEl . 

One particularly interesting class of problems con- 
cerns droplets in metallic quantum magnets. In these 
systems, the dynamics of the magnetic modes is over- 
damped because they couple to gapless fermionic ex- 
citations. In metallic (Ising) antiferromagncts, this 
dissipative environment strongly suppresses tunneling, 
and sufficiently large droplets completely freeze at low 
temperatures! 16 i 17 ' 18 The global quantum phase transi- 
tion is thus smeared f 1 ^ 

In metallic ferromagnets, the situation is further com- 
plicated because the coupling between the magnetic 
modes and the gapless fermionic degrees of freedom gen- 
erates an effective long-range spatial interaction between 
the magnetic fluctuations i 20 ' 21 i 22 This interaction which 
takes the form r~^ ld ~ v> for clean electrons and r _ ' 2d ~ 2 ^ 
for diffusive electrons, where d > 2 is the spatial dimen- 
sionality, can be viewed as a result of generic scale in- 
variance (for a recent review see Ref. H^). Understanding 
defects in nearly critical metallic quantum ferromagnets 
thus leads to the question of whether the existence and 
the properties of magnetic droplets arc influenced by the 
long-range spatial interaction. 

In this paper, we therefore develop the theory of a 
single defect coupling to the square of the order parame- 
ter in a nearly critical classical or quantum magnet with 
power-law spatial interactions of the form r~( d+CT ) with 
a > to ensure a proper thermodynamic limit. A crucial 
effect of the long-range interactions is that the tail of the 
magnetic droplet decays into the bulk region like a power- 
law of the distance as mandated by Griffiths theorem. 24 
Such a strong tail extending into the region that prefers 
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the disordered phase can be expected to be energetically 
unfavorable. To find out to what extent this hinders the 
formation of the magnetic droplet, we study the droplet 
free energy within the saddle-point approach. In the 
quantum case, we also consider the tunneling dynam- 
ics of the droplet for three cases: undamped dynamics, 
overdamped dynamics due to Ohmic dissipation, and a 
conserved overdamped dynamics as in the itinerant fer- 
romagnet. 

Our paper is organized as follows: We introduce our 
model, a classical or quantum <j> -theory with long-range 
spatial interactions, in Sec.|TTJ In Sec. lIIIl we analyse the 
free energy of a droplet within saddle-point approxima- 
tion, and we discuss fluctuations. The droplet dynam- 
ics in the quantum case is considered in Sec. HVl The 
concluding Sec. [V] is devoted to a summary as well as 
a discussion of the order parameter symmetry and the 
consequences of our results for quantum Griffiths effects. 



II. THE MODEL 

In this section we introduce our model, a d-dimensional 
Landau-Ginzburg- Wilson field theory with long-range 
power-law interactions for a scalar order parameter field 
(p. We first formulate the model for the case of a zero- 
temperature quantum phase transition, and we later dis- 
cuss the necessary changes for a classical thermal phase 
transition. The action of our quantum </> 4 -theory reads 



S — S'stat + <Sdyn> 

with the static part given by 



(1) 



/ dr J dxdyip (x, r) T (x, y) ip (y, 
| J drdxip 4 ' (x, t) . 



(2) 



Here, x and y are position vectors and r is imag- 
inary time. The bare two-point vertex, T (x, y) = 
Fni(x)(5 (x — y) + Ti(x, y), contains a non- interacting 
part and the attractive long-range interaction. The latter 
is given by 



r I (x,y) = - 7 £ 2 + |x 



(3) 



Here, 7 is the interaction strength, £0 is a microscopic 
cutoff length scale of the order of the lattice constant, 
and a controls the range of the interaction. To ensure a 
proper thermodynamic limit (an extensive free energy), 
a must be positive. Note that an additional short-range 
interaction of the usual form \Vip\ 2 can be added, if de- 
sired. As will be shown in Sec. IIII Al its contribution is 
subleading. The noninteracting part of the vertex reads 



r NI (x) = t Q + st (x) + r , 



(4) 



where to is the bulk distance from criticality^ and the 
constant Tq is chosen such that it cancels the (q = 0) 



Fourier component of the interaction (thus ensuring that 
the bulk critical point is indeed at to = 0). It takes the 
value T = Qd'ftio' 7 B(a/2, d/2)/2. Here fi^ is the surface 
of a d-dimensional unit sphere, and B(x, y) is Euler's beta 
function. St (x) is the defect potential. For definitcness 
we consider a single spherically symmetric defect at the 
origin, 



-V (|x|<a) 
(|x|>a) 



(5) 



We are interested in the case V > 0, i.e., in defects that 
favor the ordered phase. 

When discussing the quantum tunneling dynamics of 
the magnetic droplets in Sec. IIV1 we will compare three 
different dynamical actions, (i) In the undamped case, 
the dynamical action is given by 



5 £L = T E / tfq^f|£(q,Un)| 



(6) 



where 0(q, u> n ) is the Fourier transform of the order pa- 
rameter field in terms of wave number q and Matsubara 
frequency u n , T is the temperature, and c plays the role 
of a velocity of the undamped modes. 

(ii) If the magnetic modes are coupled to an ohmic 
bath, the leading term in the dynamic action takes the 
form 



S, 



(2) 
dyn 



rfq 



l<^(q, w„ 



(7) 



where a measures the strength of the dissipation, and 
there is a microscopic frequency cutoff \u> n \ < ui m i C - 

(iii) Finally, we also consider the case of overdamped 
dynamics with order parameter conservation analogous 
to the itinerant ferromagnet. The leading term in the 
dynamic action is given by 



o(3) 



(8) 



The action defined in Eqs. (JT]) to (|8]) describes a system 
close to a quantum phase transition. In order to inves- 
tigate a droplet in a system close to a classical thermal 
phase transition, we simply drop the dynamical piece of 
the action and eliminate the imaginary time-dependence 
of the order parameter field. 



III. EXISTENCE OF MAGNETIC DROPLETS 

In this section we investigate to what extent the exis- 
tence of droplets is influenced by the long-range spatial 
interaction. The basic idea is as follows: If the local po- 
tential to — V on the defect is negative, magnetic order is 
preferred on the defect even though the bulk system may 
be nonmagnetic, to > 0. Figure Q] shows a schematic of 
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FIG. 1: (Color online) Schematic local order parameter pro- 
files for defect induced droplets for short-range (a) and long- 
range (b) interactions. The dashed line depicts the defect 
potential. 



the local order parameter profile in this situation, com- 
paring short-range and long-range interactions. In the 
short-range case, the tail of the droplet profile falls off 
exponentially outside the defect i 17 ' 18 i 25 Thus the tail pro- 
vides only a subleading surface term to the droplet free 
energy. In contrast, for the long-range interaction |3]), 
the tail must take a power-law form because Griffiths 
theorem 2 ^ dictates that the magnetic correlations can- 
not decay faster than the interaction. The tail thus ex- 
tends far into the bulk region where the local potential 
is positive, and therefore leads to a large positive con- 
tribution to the droplet free energy. In this section we 
study whether this mechanism hinders the formation of 
magnetic droplets for long-range interactions. 



A. Saddle-point equation 

We start by analyzing the action |T]) within saddle- 
point approximation, focusing on the case to > (non- 
critical bulk system) because it is relevant for Grif- 
fiths phenomena. We can restrict ourselves to time- 
independent solutions because they have the lowest 
saddle-point actions (any time dependence produces an 
extra, strictly positive contribution from Sd yn )- Setting 
<p(x, r) = </>(x) and minimizing the total action with re- 
spect to this field leads to the saddle-point equation 



(to + St (x)+r o )0(x)+ M 3 (x) 



7^ (y) d y 



Co 2 



Note that the classical action discussed at the end of sec- 
tion [TT] leads to the same saddle-point equation. There- 
fore, the remainder of this section applies to both the 
classical and quantum cases. 

We have not managed to solve the nonlinear integral 
equation © in closed form. We therefore first present 
analytical results for the behavior of </>( x ) f ar away from 
the defect, and then we complement them by a numeri- 
cal solution. For sufficiently large V (such that to — V is 
sufficiently negative), we expect the order parameter in 
the droplet core, |x| < a, to be roughly constant. Grif- 
fiths' theorem 24 mandates that the droplet tail cannot 
decay faster than |x|~( d+<T '; wc therefore try the spheri- 
cally symmetric ansatz 



</,(x) 



C/|x 



d+a 



(|x| < a) 
(|x| > a) 



(10) 



with parameters 4>o and C. Note that in general, the 
ansatz (|TU)) is not continuous at |x| = a. To cure this 
unphysical behavior, there must be an intermediate re- 
gion a < |x| < a + £ m which connects the core with the 
asymptotic region in (jTUjl . We will come back to this 
point later in this section. 

We now insert the ansatz (JTOJ) into the saddle-point 
equation ^ and analyze it in the limit of large defects, 
a 3> £o, and large distance, |x| ^> a where (0 can be 
linearized in <fi. We find that the ansatz indeed solves 
the linearized saddle-point equation with the amplitude 
C given by (to leading order in a) 



_ fl d (j)oj d 

o — — a 

dt 



(11) 



(9) 



Note that C diverges when the bulk system approaches 
criticality (to — > 0) indicating that the ansatz (fT0|) is not 
valid for a defect in a critical bulk. We will come back 
to this point in the next subsection. 

To determine 0o, we now calculate the saddle-point 
action by inserting the solution (fit))) with (fTT]) into the 
action (p}. The result is the sum of a droplet core term, a 
tail term, and a core-tail interaction term. The core term 
takes the form (£ld/d)a d (j)o(to — V + u4>q/2). The contri- 
bution of the long-range interaction is exactly cancelled 
by the To-term, as must be the case for a constant order 
parameter. Interestingly, the tail term and the core-tail 
interaction term are subleading in the limit of large de- 
fects, a 3> £ - Their leading a-dependencies are a d ~ 2a 
and a d ~~ CT ~ 1 (up to possible logarithmic corrections), re- 
spectively. Finally, we have to consider the intermediate 
region a < |x| < a + £.,„ in which droplet core smoothly 
connects to the asymptotic tail. From the numerical so- 
lution of the saddle-point equation (discussed in the next 
subsection) we found that the width of the intermediate 
region is of the order of the microscopic scale, £ m ~ £o 
(at least as long as the bulk system is not too close to 
criticality; see next subsection for details). Importantly, 
£ m does not depend on the defect size a. Therefore, the 
intermediate region can at most make a surface-type con- 
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tribution to the droplet free energy, i.e., it can at most 
scale like a d_1 . 

Collecting all the terms, we find that the saddle-point 
action takes the form 

S SP = ^<t>W (to-V+ \t&) + O (a d -\ a d - 2 -) (12) 

in the limit of a large defect (a — > oo). Minimizing 5sp 
with respect to 4>o gives the optimal value 

This means, in the limit of a large defect, a droplet of 
local order starts to form as soon as the local potential 
to — V on the defect becomes negative. For finite a, the 
subleading terms in (Tl2"|) lead to a shift in the onset of 
local order that can be described by finite-size scaling in 
the usual way. 

The results f| 1 2[) and (fT"3|) are identical to the case 
of short-range interactions! 17 ' 18 > 25 We thus arrive at the 
somewhat surprising conclusion that even though the 
long-range interactions do induce a power-law tail of the 
droplet, they do not change the leading behavior of its 
free energy (in the limit of large defects), and thus do not 
hinder the existence of large droplets. 

We also note that an additional short-range interaction 
of the form \ V4>\ 2 in the static action ([2]) will not modify 
our results. Clearly, in the core region of the droplet it 
plays no role, and faraway from the core (x ^> a), it only 
produces a subleading power-law. Its contribution in the 
intermediate region can at most be of order a d_1 . 



point of the bulk rather than its critical fixed point. 
Since this stable fixed point does not have anomalous 
dimensions, the saddle-point analysis is qualitatively cor- 
rect and the decay exponent in (JTDJ) remains unchanged. 
Thus, the fluctuations only renormalize nonuniversal 
prefactors. Analogous results were found in Ref. [It] for 
the case of short-range case interaction. Note that crit- 
ical fluctuations on the defect can change the exponent 
in the relation (|T3)) close to the onset of local order at 
to — V = 0, provided the system is below its upper criti- 
cal dimension. However, this has no bearing on the form 
of the tail. 

In contrast, if the bulk system is right at the transi- 
tion, to = 0, the long-distance magnetization decay is 
controlled by the exponent rj of the critical fixed point 
via <^>(x) ~ |x|~ d+2 ~'' (because far from the defect, </>(x) 
falls off as the bulk correlation function). For a classical 
magnet with long-range interactions this fixed point was 
studied in the seminal work of Fisher, Ma, and Nickel. 26 
They found that the critical behavior is mean-field-likc 
for er < d/2 with rj = 2 — a . For a > 2 — t^sr (where rysrt, 
is the exponent of the corresponding short-range model), 
the critical behavior is identical to that of the short-range 
modeh 27 ' 28 i 29 In between, the exponents are nonclassical 
and interpolate between mean-field and short-range be- 
havior. Lets us also point out that interesting crossover 
phenomena occur when the bulk system is close but not 
exactly at the critical point. In this case the critical fixed 
point controls the magnetization decay at intermediate 
distances (of the order of the bulk correlation length) 
from the defect while the asymptotic behavior is again 
given by the saddle-point result (jTUJ) . 



B. Fluctuations 

So far, we have analyzed the magnetic droplets within 
saddle-point approximation. In this subsection we dis- 
cuss to what extent fluctuations modify the above saddle- 
point analysis. It is useful to divide the fluctuations into 
two classes, small fluctuations about the saddle-point so- 
lution and collective reorientations of the entire droplet 
in (imaginary) time. These two classes are well separated 
if the local order on the defect is properly developed, i.e., 
V — to > u. The collective reorientations determine the 
long-time quantum dynamics of the droplet. They will 
be considered in more detail in Sec. IIVI 

In contrast, small long- wavelength fluctuations could 
potentially modify the droplet profile (|10[) . in particular 
the form of the magnetization tail. To study the rel- 
evance of these fluctuations, we expand the action (|TJ) 
about the saddle-point solution and perform a tree-level 
(power-counting) renormalization group analysis. The 
results depend qualitatively on whether or not the bulk 
system is critical. 

As long as the bulk system is in its disordered phase, 
to > 0, the asymptotic long-distance decay of the droplet 
magnetization is controlled by the stable large-to fixed 



C. Numerical solutions of the saddle-point 
equation 

In this subsection we confirm and complement the 
asymptotic analysis of the saddle-point equation @ by 
a numerically exact solution. 

We study both one and three space dimensions. In the 
three-dimensional case, for a spherical defect and droplet, 
the angular integration on the r.h.s. of the saddle-point 
equation ^ can be carried out analytically leading to a 
one-dimensional integral equation in radial direction. We 
now discretize space in units of the microscopic length £o 
and fix the energy scale by setting u = 1. The resulting 
set of nonlinear equations is solved by the following pro- 
cedure: We start from an ansatz for 4> (e.g., the ansatz 
given in (|10[) ) and numerically perform the integral in the 
long-range term of ^j. We then determine an improved 
value for <f> by solving the remaining cubic equation at 
each point by standard methods. These steps are re- 
peated iteratively until the solution converges. 

In this way, we have analyzed one-dimensional systems 
with 2 x 10 4 to 2 x 10 5 points and three-dimensional 
systems with 10 4 to 10 5 points in radial direction. Wc 
have studied the cases a — 1,2,3, large defects a ^> 1 




FIG. 2: (Color online) Local order parameter ^ of a three- 
dimensional droplet as a function of distance x from the de- 
fect center for different distances to — 0.1 to 51.2 from bulk 
critica 1 " 



FIG. 4: (Color online) Local order parameter <f> of a one- 
dimensional droplet as a function of distance x from the defect 
center for different defect sizes a — 50 to 1600 (from left to 
right). 
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FIG. 3: (Color online) Amplitude C of the asymptotic power- 
law decay of the droplet tail for the system shown in Fig. [2] 
The solid line is the theoretical prediction, Eqs. (|ll[l and (|13|) . 

and various values of to, V and 7. For weak long-range 
interactions and away from bulk criticality, our procedure 
converges rapidly. With increasing 7 and decreasing to, 
the convergence becomes slower. However, in all cases, 
our self-consistency cycle eventually converges, giving us 
a numerically exact solution of the saddle-point equation. 

We now present and discuss a few characteristic results 
from these calculations. In Fig. [2J we show saddle-point 
solutions for d = 3, a = 1 and different values of the 
distance to from bulk criticality. In agreement with the 
analytical predictions of the last subsection, the order 
parameter is essentially constant on the defect. For large 
|x|, the droplet tail falls off with the predicted power-law 
cj) ~ |x|~( d+<T) = |x|~ 4 for all values of t . The amplitude 
C of this power-law decay is analyzed in Fig. [3) As pre- 
dicted in the last subsection, for small to, C behaves like 
1/to (the small deviations are the lowest to stem from 
the fact that in these cases, 10 5 sites is not sufficient to 
reach the asymptotic regime). 

Figure Q] shows the dependence of the droplet profile 



FIG. 5: (Color online) Amplitude C of the asymptotic power- 
law decay of the droplet tail for the system shown in Fig. [4] 
The solid line is the theoretical prediction, Eqs. (111[) and (|13|) . 



on the size a of the defect for a system with d = a = 1 . 
For all a, the asymptotic decay of the droplet tail takes 
the predicted power-law form, <j> ~ |x| — C d +°") = |x|~ 2 . 
This figure also shows that the width £ m of the interme- 
diate x-region which connects the droplet core with the 
power-law tail does not change with a as discussed in the 
last subsection. (This becomes even more obvious when a 
linear rather than the logarithmic x-scale is used.) More- 
over, in agreement with (|13p. 4>o does not depend on a. 
The amplitude C of this power-law decay is analyzed in 
Fig. [SI In agreement with the theoretical prediction pip , 
the amplitude grows linearly with the defect size a. 

We have performed analogous calculations for other 
parameter sets, including varying to for a — 
300, 600, 1600 as well as varying a for V — 30. In 
all cases, we have found excellent agreement with the 
predictions of the asymptotic analysis of Sec. IIII CI 
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IV. TUNNELING DYNAMICS 

In this section, we study the tunneling dynamics of 
a single droplet at a zero-temperature quantum phase 
transition. Our approach starts from the pioneering work 
of Callan and Coleman 3 -^ and Leggett and coworker a 31 i 32 
(in the case of dissipative dynamics). In the following 
subsections we separately discuss the droplet dynamics 
for the three dynamical actions given in Eqs. © to ©, 
starting with the undamped case. 



A. Undamped magnet 

Following Callan and Coleman^, the tunneling rate 
between the up and down states of the droplet (i.e., 
the tunnel splitting of the ground state) can be esti- 
mated from the action of instanton-like saddle-point solu- 
tions y(x, t) fulfilling the boundary conditions <p(x, r) —> 
±</>(x) for t — > ±00. In principle, several processes con- 
tribute to the overall tunneling rate. In the simplest 
one, the droplet retains its shape while collapsing and 
reforming. 33,34 A competing process consists of the nu- 
cleation of a domain wall that then sweeps the droplet^ 

We start by considering the collapse-and-reformation 
process which can be described by an ansatz 



<j> (x) r) (t) 



(14) 



with </>(x) being the static saddle-point solution of section 
Mil and T){t) —> ±1 for t — » ±00. Inserting this ansatz 
into the action (TTJ) and integrating over the spatial vari- 
ables yields the following excess effective action (above 
the time- independent solution 77 = 1) 



AS™ 



n 



U 2 a d / dr 



4>pU 
2 



1 



1 



V 



dr] 

>, 

(15) 

to leading order in the defect size a. The saddle-point 
instanton solution of this action can be found exactly. It 
takes the form rj (r) = tanh (r/ro), with t^ 2 = c 2 (/)qu/2. 
The resulting instanton action reads 



(i) 



AS; 



giving a tunnel splitting 



(1) 



40o 

3d 



LO e 



u4>oa d T 



-AS K 



(16) 



(17) 



The "attempt frequency" uq can be determined by stan- 
dard quantum tunneling considerations^ 2 - from the fluc- 
tuations about the instanton solution. Importantly, the 
tunneling rate decays exponentially with the volume of 
the droplet. Equations ([15]) to (TIT)) are in complete agree- 
ment with the corresponding results for the case of short- 
range interactions ) 17 ^ 8 ' 36 reflecting the fact that the lead- 
ing terms of the instanton action stem from the droplet 
core rather than the tail. 



To discuss the contribution of the moving domain wall 
processes to the tunneling rate, we use the ansatz 



0(x)r/(r — x/v) 



(18) 



which describes a domain wall that sweeps the droplet in 
x-direction with velocity v. r\ describes the domain wall 
shape and fulfills the boundary conditions rj(z) — > ±1 for 
z — > ±00 as before. Inserting this into the action ((T|) 
gives the same effective action (fT5"j) plus one additional 
positive term from the spatial dependence of rj (this 
term corresponds to the domain wall energy). There- 
fore, the minimal action for this process is bounded by 
(I16p . and to exponential accuracy the corresponding tun- 
nelling rate cannot be larger than (|17p . This is in qualita- 
tive agreement with earlier results for short-range inter- 
actions. Stauffer— estimated the tunneling rate of a do- 
main wall within quasiclassical WKB approximation and 
found that it depends exponentially on the droplet vol- 
ume. Senthil and Sachdev 3 -^ estimated the tunnel split- 
ting of a locally ordered island in a transverse-field Ising 
model using perturbation arguments. Again, the result 
(which should contain all possible processes) is exponen- 
tially small in the droplet volume. 



B. Overdamped dynamics 



We now consider overdamped dynamics with the ac- 
tion S^ yn as given in ([7]) . Inserting the ansatz <p (x, r) = 

4> (x) rj (r) into S\' and integrating over the spatial vari- 
ables gives the following contribution to the effective ac- 
tion 



drdr 



, a (v(t) - n{ T ')f 



2tt 



(t- 



(19) 



to leading order in the defect size a. The other terms are 
as in eq. (|15[) . A straight forward saddle-point instan- 
ton analysis of the effective action analogous to the last 
subsection fails because the interaction of the trajectory 
7](t) at large positive times with the trajectory at large 
negative times causes a logarithmic divergence. Follow- 
ing Refs.HH and[H, the calculation therefore proceeds in 
two stages. 

In the first stage, we introduce a low-frequency cutoff 
lu c in the dynamic action ([7]). This changes the interac- 
tion kernel in (TT91). 



1 



1 + 2lu c \t - t'\ 



{t-t'Y {t - t'Y {1 + u c \t 



(20) 



and removes the divergence. We have not been able to 
solve analytically for the instanton but we have used the 
ansatz T)(t) — tanh(r/ro) with variational parameter to. 
Minimizing the effective action AS^ + AS^ with re- 
spect to tq gives 



3a 



1 



2n^u 
9c 2 a 2 



(21) 
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In the limit of weak dissipation, a — > 0, we recover the 
result for undamped dynamics while strong dissipation, 
a — > oo, leads to To = 6a/ '(ircj^u). The resulting instan- 
ton action can be expressed in terms of the dimensionless 
dissipation strength paramete r 32 ' 38 



—(j) z a d a 
ira 



(22) 



We note that a is proportional to the defect volume 
a d . (Analogous results have been obtained for dissipa- 
tive random quantum Ising models i 39 ' 40 ) The dissipative 
part of the instanton action reads 



AS-, 



(2) 



-aln(w c ) + /(a) 



(23) 



where the function f(a) is given by f(a) — ca + 0(a 2 ) 
for weak dissipation and /(a) = — alna + c'a + 0(ln(a)) 
for strong dissipation, c and c' are constants of order 
one. For comparison we have also studied a piecewise 
linear ansatz for tj(t). The resulting instanton action is 
identical to ([23]) except for different numerical values of 
the constants c, c'. At the end of the first stage of the 
calculation we thus obtain the bare tunnel splitting 



(2) 
bare 



UJ e 



-AS} 



(24) 



In the second stage of the calculation the resulting dis- 
sipative two-level system is treated using renormaliza- 
tion group methods.— It is well-known that instanton- 
instanton interactions renormalize the tunnel splitting, 
yielding 



,(2) 



J {2) 
bare 



,(2) 



*/(!-«) 



(25) 



This also eliminates the unphysical dependence of the 
tunnel splitting on the arbitrary cutoff parameter u> c . 
We thus find that the smaller defects with a < 1 con- 
tinue to tunnel, albeit with a strongly reduced rate. The 
larger defects with a > 1 cease to tunnel, i.e., they are on 
the localized side of the Kosterlitz-Thouless phase tran- 
sition of the dissipative two-level system. These results 
are in qualitative agreement with the case of short-range 



interactions 



17.18 



C. Conserved overdamped dynamics 

Finally, we consider the case of overdamped dynamics 
with a conserved order parameter as given by the dy- 
namic action ([8]). Such an action arises, e.g., in the case 
of an itinerant quantum ferromagnet. 

Order parameter conservation requires some care in 
discussing the dynamics of our locally ordered droplet. In 
particular, the homogeneous magnetization J dx(/?(x, r) 
must not be time dependent. Therefore, the product 
form (f (x, t) — 4> (x) -q (r) with </>(x) the static solution 
of section IIIII is not a suitable ansatz in this case. This 



can be fixed (in a crude way) by subtracting a constant 
from the droplet profile, <fi'(x) — </>(x) — const such that 
the q = Fourier component is cancelled. The ansatz 
tp (x, r) = </>' (x) r\ (r) then provides a variational upper 
bound for the instanton action. 

Inserting this ansatz into ([5]) and carrying out the inte- 
gral over the spatial variables leads to a dissipative term 
in the effective t](t) action with the same functional form 
as ([T9|) . The prefactor and the resulting dimensionless 
dissipation strength a, however, are different. To leading 
order in the defect size a, we find 



32^a 3 a c /(37r) 



(d = 3) 
(d = 2) 



(26) 



In general dimension d > 2, the dimensionless dissipation 
strength is now proportional to a d+1 instead a d . The ex- 
tra factor a compared to the non-conserved case in Sec. 
HVBI can be understood as follows. To invert the magne- 
tization of a droplet of linear size a, magnetization must 
be transported over a distance that is at least of order a 
(because the order parameter conservation prevents sim- 
ply flipping the sign of the magnetization on the defect) . 
This involves modes with wave vectors of the order of 
q ~ 1/a. Since the dissipation strength in ([5]) is inversely 
proportional to q, we expect an additional factor a in the 
effective action. This argument strongly suggests that 
this extra factor is not an artefact of our simple ansatz 
for if (x, r) but correctly reflects the physics of the con- 
served order parameter case. 

In all other respects, the calculation proceeds as in the 
non-conserved case in Sec. I IVES I The resulting dynamic 
behavior of the droplets depends on the value of the 
dimensionless dissipation strength parameter a. Small 
droplets (a < 1) still tunnel while the larger ones (a > 1) 
freeze. Because a is now proportional to a d+1 the tun- 
neling of large droplets is even more strongly suppressed 
than in the non-conserved case. 



V. DISCUSSION AND CONCLUSIONS 

To summarize, we have studied the physics of a single 
defect coupling to the square of the order parameter in a 
nearly critical system with long-range spatial interactions 
of the form r -( d + a ) with a > 0. Such a defect can in- 
duce the nuclcation of a magnetic droplet while the bulk 
system is still in the nonmagnetic phase. Due to the long- 
range interactions, the droplet magnetization develops a 
long power-law tail, i.e., at large distances r from the de- 
fect, it decays like r -( d + cr ) i n agreement with Griffiths' 
theorem<2i Nonetheless, the droplet free energy is domi- 
nated by the core (on-defect) contribution while the tail 
contribution is subleading in the limit of large defects. 
Therefore, droplets will nucleate on large defects as soon 
as the local potential (the local distance from criticality) 
becomes negative, in complete agreement with the case 
of short-range interactions. Our explicit calculations of 
the droplet magnetization profile have been performed 
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within saddle-point approximation, but as long as the 
bulk system is noncritical, fluctuations do not change the 
functional form of the droplet. They only renormalize 
nonuniversal parameters. 

In addition to the existence of the magnetic droplets, 
we have also investigated their dynamics. As is well 
known^ in the case of a classical (thermal) phase transi- 
tion, the droplet cannot order statically. Instead, it fluc- 
tuates between 'up' and 'down' due to thermal fluctua- 
tions. For a zero-temperature quantum phase transition, 
the behavior is potentially different, depending on the 
form of the dynamic action. We have studied three cases. 
In the absence of dissipation, even very large droplets can 
always tunnel, but with a rate that decreases exponen- 
tially with the droplet volume. This changes in the pres- 
ence of (Ohmic) dissipation. The qualitative behavior 
now depends on the dimensionless dissipation strength 
a. For a < 1, the droplet still tunnels albeit with a 
further reduced rate while for a > 1, tunneling ceases 
and the droplet magnetization becomes static. For over- 
damped dynamics without order parameter conservation, 
a is proportional to the volume of the droplet core. Thus, 
sufficiently large droplets always freeze in agreement with 
Refs.[lllllll. 

In the case of overdamped dynamics with 
order parameter conservation as in the itinerant quan- 
tum ferromagnet, the dissipation effects are further en- 
hanced because the dimensionless dissipation strength a 
for a droplet of linear core size a is proportional to a d+1 
rather than a d . 

Let us comment on the order parameter symmetry. 
Our explicit results have been for the case of a scalar 
(Ising) order parameter. However, the analysis of the 
droplet existence in Sec. |TTT] relied on saddle-point ar- 
guments and thus applies equally to continuous O(N) 
order parameters with TV > 1. In contrast, to generalize 
the discussion of the dynamics in Sec. |IV] to such order 
parameters, other types of fluctuations (rotational ones) 



must be considered. 

We also emphasize that we have discussed the case of 
an isotropic attractive long-range interaction. Droplet 
formation dominated by oscillating and/or anisotropic 
interactions such as the dipolar or the RKKY interactions 
is likely of different type and not considered here. 

Finally, we briefly discuss the consequences of our re- 
sults for the (quantum) Griffiths effects in systems with 
long-range spatial interactions. Because the power-law 
magnetization tail only makes a subleading contribution 
to the free energy of a magnetic droplet, such droplets 
can form on rare (strongly coupled) spatial regions of 
the disordered system essentially in the same way as in 
the case of short-range interactions. Therefore, as long 
as droplet-droplet coupling can be neglected, the Grif- 
fiths effects should be identical to those in short-range 
interacting systems. However, it is clear that the droplet- 
droplet coupling is more important for long-range inter- 
actions than for short-range ones. This means, it must be 
considered for lower droplet density and, in the quantum 
case, for higher temperatures. The complicated physics 
caused by the coupling of several droplets is beyond the 
scope of this paper. Recently, it has been argued 4 ^ that 
this coupling can qualitatively change the Griffiths effects 
at least in some cases. A complete understanding of this 
phenomenon remains a task for the future. 
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